Lieb-Robinson bounds for commutator-bounded operators 
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The principle of locality is at the heart of the foun- 
dations of all modern physics. In quantum field the- 
ory, the principle of locality is enforced by an exact light 
cone. Whenever two (bosonic) observables are space- 
like separated, they have to commute, so that neither 
can have any causal influence on the other. In ordi- 
nary quantum mechanics, no explicit request for local- 
ity is imposed, and it is, in principle, possible to sig- 
nal between arbitrarily far apart points in an arbitrarily 
short time. Nevertheless, a simple perturbation analy- 
sis shows that such an influence must decay exponen- 
tially with the distance between the observables. The 
seminal work by Lieb and Robinson [lj has made this 
statement rigourous for nonrelativistic spin systems. 
In essence, it states that any quantum system whose 
Hilbert space is composed of a tensor product of lo- 
cal, finite-dimensional Hilbert spaces and whose Hamil- 
tonian is the sum of local operators will have an ap- 
proximately maximum speed of signals. Here, local just 
means that every operator has as a support the tensor 
product of few degrees of freedom. The approximation 
consists of the fact that outside the effective light cone 
there is an exponentially decaying tail. 

Recently, Lieb-Robinson bounds (LRBs) have received 
renewed interest in both the fields of theoretical con- 
densed matter and quantum information theory |2]- 
[l0l[T3WT7l . In particular, the LRB has been used to prove 
that a nonvanishing spectral gap implies an exponential 
clustering in the ground state |6 , 8 13J. Further develop- 
ments can be found in |9|, where the LRB is used also to 
argue about the existence of dynamics. The LRB has also 
been instrumental in the recent extension of the Lieb- 
Schultz-Mattis theorem to higher dimensions ITPTI [I4l . 
In Il2l , it has been shown how the Lieb-Robinson 
bounds can be exploited to find general scaling laws 
for entanglement. In [2J, these techniques have been ex- 
ploited to characterise the creation of topological order. 
The locality of dynamics has important consequences on 
the simulability of quantum spin systems. In 1201 I2T1 it 
has been shown that one-dimensional gapped spin sys- 
tems can be efficiently simulated. A review of some of 
the most relevant aspects of the locality of dynamics for 
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quantum spins systems can be found in |15|. Other de- 
velopments of significant interest include 11611171 which 
show that it is possible to entangle macroscopically sep- 
arated nanoelectromechanical oscillators of the oscilla- 
tor chain and that the resulting entanglement is robust 
to decoherence. Such a system is of great interest for its 
possible application as a quantum channel and as a tool 
to investigate the boundary between the classical and 
quantum worlds. 

The LRBs have found a more exotic use in the field 
of emergent gravity, where one wants to study local- 
ity, geometry and Lorentz symmetry as emergent phe- 
nomena Il22l l23l . An example of the usefulness of the 
Lieb-Robinson bounds can be found in [24J, where it 
was shown that in spin systems with emergent electro- 
magnetism 1 19 1, the speed of light is also the maximum 
speed of signals, without imposing from the beginning 
any Lorentz invariance. This raises the issue whether 
even Lorentz invariance could be emergent. 

One problem with the Lieb-Robinson bounds is that 
it is difficult to obtain bounds for unbounded Hamilto- 
nians. In the usual Lieb-Robinson settings, the Hamil- 
tonian must be a sum of local bounded operators. If the 
unbounded terms in the Hamiltonian are completely lo- 
cal, that is, if they are on-site terms, it is possible to 
prove a Lieb-Robinson theorem using the usual tech- 
nique 1 10 1 . In the specific case of coupled harmonic os- 
cillators on a graph with local interactions, it was proven 
in J3J that the Lieb-Robinson bound is valid for canon- 
ical and Weyl operators and a proof for a generalisa- 
tion to general operators is outlined. Algebraic suppres- 
sion (instead of the usual exponential suppression of the 
Lieb-Robinson Bound) is shown to result from nonlo- 
cal algebraic interactions. As an interesting corollary, 
1 3 1 shows how the approximate locality implied by the 
Lieb-Robinson bound becomes exact in the continuous 
limit for the Klein-Gordon field. 

In this paper, we will show how one can find a bound 
to the maximum speed of interactions in the case of a 
class of unbounded spin Hamiltonians. It is not true 
that for any unbounded Hamiltonian, a Lieb-Robinson 
bound exists (181 . Here, we want to show that one can 
derive a Lieb-Robinson bound if the Hamiltonian is the 
sum of local operators, whose commutators are bounded. 
Therefore, there is no necessity for even the nonlocal 
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terms to be bounded, as long as their commutators 
are. More specifically, we show that for quantum sys- 
tems whose Hilbert space is the tensor product of local 
Hilbert spaces associated with vertices and edges of a 
graph, if the Hamiltonian is the sum of local operators 
<&i, each with a support on a region of the graph with a 
diameter less than a fixed number i?,if each of these local 
operators $i is noncommuting with less than v other lo- 
cal operator terms of the Hamiltonian <£>.,, and if for any 
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A (v LR t - d) 
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where M is a constant and vlr, the limit on the speed of 
propagation of information, depends only on local op- 
erators of the Hamiltonian as it is they who affect the 
propagation. This bound can be generalised to any local 
observables Op and Oq with supports P and Q respec- 
tively, that satisfy the following local observable opera- 
tor conditions: (i) The graph distance d separating P and 
Q is greater than R. (ii) The number of terms of the 
Hamiltonian whose support has nonempty intersection 
with P is np < oo. (Hi) There exists Fp and Fq such 
that for all terms <£>i and <i>j of the Hamiltonian the in- 
equalities || [Op, $(]|| < F P K, WpQ^iW < F Q K and 
|| [Oq, [«&,-, || < FqQ are satisfied. The generalised 
bound is then 



\\[Op(t),0< 

< MFpF Q n P {n P + 1) exp [a (v LR t - d) ] . (2) 

To motivate our discussion, let us start by the most 
trivial example: Consider the case of a Hamiltonian H = 
^2 hi which is composed of a sum of local terms hi which 
are commuting, such as the quantum Ising model with- 
out the transverse field. In such a case, there is simply 
no propagation of signals: Indeed, for any local operator 
A we have A {t) = e ltH A e~ iHt = e itBA A <r iHAt , 
where H A = J2i-[hi o A ]^o ^" smce there is a finite num- 
ber of hi in H A , and they are of finite range; OA{t) is 
also strictly local for arbitrary long times t, irrespective 
of the norm of the hi operators. This suggests that it is 
desirable to find Lieb-Robinson bounds in terms of the 
norm of the commutators rather than the norm of the 
local terms hi. 

Let us outline a simple example. Consider a system of 
parallel quantum wires. We place fermions on the wires, 
and these are usually described by one-dimensional 
Luttinger liquids, and have approximately a linear dis- 
persion relation. We place a density-density interaction 
between the wires. Labeling the wires by the index j, the 



system can be described by the following Hamiltonian: 
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which is commutator bounded in the sense of this paper 
as long as both and \d Xi d x .V\ are bounded. An- 

other example involves a generalised Dicke model, de- 
scribing an array of spins interacting with a boson field 



two of these operators we have || < K and for Vla 

any three operators || [$_,-, $fc]]|| < Q for two positive 
numbers K and Q, then we have that for any two local 
operators <£>i and $i which are terms in the Hamiltonian 
and whose support is separated by a graph distance d, 
that 
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ib„+i) (4) 



where b n are boson creation operators and of is the nth 
spin. It is easy to check that in this case the commu- 
tator [h n ,h n+ i] = — 2iofof +1 is bounded. [In fact, this 
particular Hamiltonian can also be written as a sum of 
commuting terms h n — b n (a^ — i(Tn-i) + h.c.]. 

We consider Hamiltonians that are the sum of two dif- 
ferent types of operators $ and $i : 



H = 



ho% 



E h M- 



(5) 



Here, So, Si are two sets of labels, ho and hi are two 
coupling constants, and [$j,$o] = = for 

every As an example, consider the Ising model. 
Then $q = crfaf +1 and $\ = of. We call the subgraph 
which is the support of the operator $™, T(q, m) and for 
(a, b) £ {0, l} 2 , we define 
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We consider what we will refer to as commutator- 
bounded i?-local quantum systems. For such systems, 
the commutators and the commutators of commutators 
of operators of the Hamiltonian are uniformly bounded 
while the operators themselves may be unbounded and 
the operators of the Hamiltonian have support on sub- 
graphs of size less than R, for R an arbitrary natural 
number. Explicitly, the diameter of all T(q, m) is less 
than R and for any three operators <& % a , and ap- 
pearing in the Hamiltonian with the coupling constants 
h a , h b , and h c , we have that /i n /i6||[$*,$^]|| < K and 

/lafrb/icllll^ii^bl^clll < Q> where K and Q are positive 
real numbers. Note that a bounded system, which is 
uniformly bounded by K, must satisfy K < 2K 2 as well 
as Q < 4\K 3 , and thus boundedness implies commuta- 
torboundedness . 

By taking the derivative of Eq.(|6} with respect to t, 
we obtain {tQ b j {t))' = [[-iH(t),¥ a (t)}, after keep- 
ing only the terms in H(t) which do not commute with 
&a(t), and after some algebra, we get [here and in the 
following by a + 1 we mean a + 1 mod(2)] 



{K\\t))' = [K\\t)\ 



+(-ih a+ i) 



E 



ih a +i ^ 
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where, if i G So, then Zi is the finite subset of Si such Taking the second derivative, and using the fact that 
that j G Zi r(0, i) f] r(l, j) 7^ and vice versa for <i>^(t)] = 0, we obtain, after some algebraic ma- 

i G Si. nipulation, 



iaeZij isez, 2 jae^ i 2 eZi 1 

-^+1 ^ ^ [[^W.^iWU^jW,^]]. (8) 



Defining the following unitary operator U2+z{t) = 
e V 1 2 / and its associ- 



ated unitary evolution T 2+3 (t)0 = U 2+3 (t^ OU 2+ z{t), 
integrating Eq.J8|, and taking the norm, we obtain, af- 
ter some manipulations, 



\K\ j '(t)\\<h a+ i 53 ii[[^,^ 2 +1 ],^] 



fdsLhi+iK 53 53 



2h a h a+ iK 53 53 hk 3 ( S ),$ 



(9) 



where, we used the fact that || [$ a (i), ^a+i(*)lll — ^ • By integrating Eq.||9|, we get 



\K^(t)\\<\\[K\H]\\+ha + l E \ma,K 2 + llH}\\t 

i 2 £Z n 



f ds f di [2hi + iK 53 53 11^^(011+2^0+1*- E E wk 

J ° J ° \ i2£Z zl j 2 £Z il i2£Z il i 3 eZ i2 



(10) 



Since the commutators are bounded, we have H^l 1 ,^]!! < A and 3>o+i1j - <3 for some A', Q > 0. Noting 

that || , || = if r(a, ii) and r(6, j) do not overlap and || [[$£ , $^ 2 +1 ] , $j]|| = if j) does not overlap with 
either T(a, ii) or T(a + 1, i^), we see that Eq.{lO) implies 



\K b i(t)\\<KSi+ 53 h a+ iQ5{ iUi2 t 
i 2 ez il 



where we have used the following symbol: 



fds f di\2hi + i 53 53 11^(011 + 2^+! 53 E uncoil 



(11) 



sk = j i if r( fl „i)n%,fc) ^0, 

otherwise. 



(12) 
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Solving for j | K% {t)\\, we find 

\\K b j (t)\\<KSi+ ]T Qha+iSi lUi2 t+ E E 2h a+iK8 3 j2 ^+ £ E 2^ tt+1 ^i| 

+ E E E 2^ +1 ^^ 2Ul3 |+ E E E 2w£*g*u|| 

i2&Zi 1 i 3 £Zj 2 i 2 £Zi 1 i 3 £Zi 2 i4GZi 3 



+ f dv f du J u d S f di((2h a+l h a fK 2 E E E uncoil 

\ »2eZ il i 3 £Z j2 j 3 ez j2 

+ (2h 2 a+1 )(2h a+1 h a )K 2 E E E E 

l2£Zi 1 j2&Zi 1 l3^Zj 2 l4^Zi 4 

+ {2h a+1 h a fK^ E E E ik+voii + (2hD(2h a+1 h a )K 2 E E E W(0ii 

t2^Zi 1 t2£Zi 1 is£Zi 2 i^Zi 3 %2^Zi x i^GZi 2 %2,^Z\ 2 3z^Zi 2 

(13) 



Iterating this procedure we obtain by induction 

\\Kj{t)\\ < y^M^"LJ_ c „, (i4) 

ra=0 

where M = \/2 max{^, ^} x max{ |, 1} x max{^, 1} 
and where c„ is a combinatorial factor counting the 
number of linking operator chains of n operators be- 
tween T(a, ii) and T(b, j). What we call an operator 
chain is heuristically a sequence of intersecting opera- 
tors linking the initial and final operators. The process 
of constructing the sequence of operators forming the 
chain is as follows: The 2j th operator in the chain has to 
be noncommuting with the (2j — l) th one. This imposes 
that the two consecutive operators of a chain have to (1) 
be of a different interaction type and (2) have overlap- 
ping support. For the odd-numbered operators, there is 
an extra choice: The (2k + l) th operator can be an op- 
erator that does not commute with the 2k th operator (as 
for the even case) or an operator that does not commute 
with the (2k — l) th operator. That is, even operators in 
the sequence must be noncommuting with the previous 
operator in the sequence and odd operators in the se- 
quence must be noncommuting with either of the two 
previous operators in the sequence. From the recursive 
Eq.(|TT), we see that if we start with an operator ii of 
type a, the next operator in the chain must be an opera- 
tor i 2 £ Z ix of the other type (a + 1). The fact that it is 
in Zi 1 means that its support overlaps with it's, which 
is similar to what was found in the bounded case. How- 
ever if we look at the operator that comes after ii, we see 
that we have two distinct possibilities. The first [second 
double sum under the integrals of Eq.lfTTjl] is that it can 
be i3 € Zi 2 an operator of type a (different from a + 1) 
whose support overlaps with ^'s; if this were the only 
possibility, we would have exactly the same situation as 
we had for bounded systems, but the first double sum 



under the integrals of Eq.|[iT) adds another possibility. 
That second possibility (first double sum) is choosing an 
operator j% G Zi x after the operator i 2 , 32 is, like 12, an 
operator of type a + 1 which (by virtue of being in Z i± ) 
has a support overlapping it's. To find the next opera- 
tor after that, we reiterate Eq.(|TT) and thus, like the first 
operator after it, we need to choose an operator which 
is of a different type than the last one (be it j% or i 3 ) 
and has overlapping support with the last one; thus, at 
this point, we cannot "change our mind". We can thus 
see the process of building the chain as, for every two 
choices, we must choose an operator that links with the 
previous one, but every other choice, we can also choose 
an operator that links to the penultimate one instead. 

Because every two choices in building up the chain 
we must choose an operator of a different type than the 
previous one, in the end, the chain contains the same 
number of operators of type as of type 1 (plus or minus 
one). This means that there will be the same number 
of factors of ht as of ho in every term; hence, we can 
pull them out of the sums over chains and simply write 
an overall factor of y/hohi in front while passing from 
Eq.g3[ to Eq.gjj). 

Furthermore, we can always find a bound of the fol- 
lowing type for c n : 

c n < M 7 V ( t- d \ (15) 

where A is an arbitrary positive real number. This is be- 
cause the r(a, i)'s have a diameter of R or less. Hence, 
if the distance d between the initial and final points is 
greater than R n , then there are no possible linking oper- 
ator chains of n local operators between the initial and 
final points. Furthermore, since at every odd step along 
the chain there is a choice of at most v local operators 
to choose from for the next operator in the chain and 
at every even step there is at maximum 2v operators to 
choose from, there is, at most, (^/2v) n possible local op- 
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erator chains of n operators starting from any given po- 
sition. Thus, we certainly have that 

c„ < VTiy n e x(Rn ~ d \ (16) 

where A is arbitrary. Using Eq.|[T5| with Eq.(|l4}, we ob- 
tain the LRB of Eq.Q, 

IPl (*),<*#)] II < tiexpXfahohiK^eh - d), (17) 

where M — MM. To obtain the generalisation to local 
operators O p and Oq satisfying the local observable op- 
erator conditions enounced in the introduction, we in- 
troduce K^{t) = [$*(£), O Q (0)}. Using exactly the same 
procedure used to obtain Eq.|[9j, we get 

\\[O P (t),O Q (0)]'(t)\\ < 

+ 2h a h a+1 KF P £ E W R a( s )\\) 

<2ma,x{h 2 Q ,h 2 1 }np(n P + 1) f ds\\K%(s)\\, (18) 

Jo 

where Zp is the set of of labels of the terms of 
the Hamiltonian which do not commute with Op 

(\\Z P \\ = n P ), where k is such that f Q ds\\K%(s)\\ = 

ma,x ie z P Jo rf s ll-^'a( s )|| an d where, unlike in Eq.j^, the 
terms containing no integrals do not appear here be- 
cause of the condition that d > R. \\K%(s)\\ can then be 



treated in exactly the same way as \\K^ (s)|| was, with 
the only exception that while bounding the final com- 
mutators [i.e., when we place the 8 of Eq.|[T2|], we will 
need an extra factor of Fq. Thus, we obtaiivEq.|j2]|: 

\\[O P (t),O Q (0)}\\ < 

F P F Q n P (n P +l)M exp X (2\/h hiKjeh - d).(19) 

Optimising for A, we have that the Lieb-Robinson speed 
is thus 

v LR = 2^h h 1 K. (20) 

We can compare Eq.|[T9| with the bound obtained for 
Hamiltonians composedof bounded local operators and 
for bounded local observables Op and Oq, which is EUl 

\\[O P (t),O Q (0)}\\ < 

\\Op\\\\0 Q \\n P M exp X - <l). (21) 

To summarize, in this paper, we have shown that a 
Lieb-Robinson bound exists for those Hamiltonians that 
are the sum of local operators whose commutator is 
bounded. This allows for treating a class of systems with 
unbounded operators. 
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